NMR study of small molecule adsorption in MOF-74-Mg 
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We calculate the carbon nuclear magnetic resonance (NMR) shielding for CO2 and the hydrogen shieldings 
for both H2 and H2O inside the metal organic framework MOF-74-Mg. Our ah initio calculations are at the 
density functional theory level using the van der Waals including density functional vdW-DF. The shieldings 
are obtained while placing the small molecules throughout the structure, including the calculated adsorption 
site for various loading scenarios. We then explore relationships between loading, rotational and positional 
characteristics, and the NMR shieldings for each adsorbate. Our NMR calculations show a change in the 
shielding depending on adsorbate, position, and loading in a range that is experimentally observable. We 
further provide a simple model for the energy and the NMR shieldings throughout the cavity of the MOF. 
By providing this mapping of shielding to position and loading for these adsorbates, we argue that NMR 
probes could be used to provide additional information about the position at which these small molecules 
bind within the MOF, as well as the loading of the adsorbed molecule. 



I. INTRODUCTION 

Metal Organic Frameworks (MOFs)P^ have become 
very popular over the last decade, as is evident by 
their prevalence in recent studies and generous review 
in the literature.'^K^ This interest is largely due to the 
wide range of applications that have been identified for 
M OFs, ranging from molecular gas storage (C^a EMMI 



N2,[ia^C02,™I^Hj^ia) to gas separation,ffil&'drug 
delivery^^ sensing, ^^ ^^ catalysis ,^^ and photoca- 
talysis.l^^^ The utility of MOFs comes from their in- 
teractions with small molecules such as H2, CO2, and 
H2O. It is thus critical to understand the details of the 
binding process when a small molecule is adsorbed into 
the MOF. To this end, IR and Raman spectroscopy have 
been used extensively to study small molecule adsorption 
in MOFs,E2^S5] ]3^t it can be difficult to determine where 
the reactive sites reside under different loading scenarios. 
Also, for these probes the strong signals originating from 
the vibrational modes of the gas present in the experi- 
ment chamber and the MOF itself can often dominate 
the spectrum, making the analysis of the weak adsor- 
bate signal challenging. In the following, we argue that 
NMR — which has already been used successfully to study 
MOFs in a number of case^^^^^ — can be used to facilitate 
a more detailed understanding of the static behavior of 
MOF/adsorbate interactions and binding under various 
conditions. In particular, we show that NMR can pro- 
vide information about the position at which these small 
molecules bind within the MOF, as well as the loading of 
the adsorbed molecule. 

In this work, we consider the particular MOF struc- 
ture MOF-74-Mg,^^ which has been shown to have very 
high efficiency when capturing C02r^a key property for 
gas separation and storage applications. For the small 
molecule adsorbed in the MOF we consider H2, CO2, 
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and H2O. The first is obviously interesting for hydrogen- 
storage applications, while the second one is of interest in 
carbon-capture applications. However, water by itself is 
not necessarily interesting for applications, were it not for 
the fact that it strongly impedes the performance in the 
first two cases. In other words, the presence of water, due 
to its strong binding characteristics, decreases the per- 
formance of MOFs in hydrogen-storage and CO2 capture 
applications, such that its careful study is warranted.^ 
The three molecules investigated in this study, i. e. 
H2, CO2, and H2O, bind inside the MOF through phy- 
sisorption. Thus, it is apparent that the proper inclu- 
sion of van der Waals interactions is crucial for the en- 
tire study. Therefore, we use density functional the- 
ory (DFT), utilizing the van der Waals including func- 
tional vdW-DF ^^ ^^ to map the shielding of an adsorbed 
molecule within MOF-74-Mg to various characteristics. 
This truly non-local exchange-correlation functional has 
already successfully been applied to study s mall molecule 
adsorption in a variety of MOFs.E^EIHSSlHlII 



II. COMPUTATIONAL DETAILS 

The interaction of II2, CO2, and H2O with MOF- 
74-Mg was studied using DFT with vdW-DF as imple- 
mented in QuantumEspresso.^ It is weh known that 
binding distances will usually be slightly overestimated 
using vdW-DF.^^ We used the primitive rhombohedral 
unit cell of MOF-74-Mg with space group R3 and 54 
atoms. The initial geometry of MOF-74-Mg was relaxed, 
fixing the lattice parameters according to the experi- 
mental values of a = 15.117 A and a = 117.742°. ^^^s A 
complete volume relaxation for all loadings considered 
in this paper would have been extremely computation- 
ally expensive, so that we fixed the lattice constants 
to the experimentally measured ones. In our testing, 
for the expected worse case of 12 II2O molecules ad- 
sorbed in MOF-74-Mg, we find that there is only a 1% 



change in the channel diameter when a fuh relaxation 
is performed. For each adsorption case we have relaxed 
the internal coordinates until the total force was below 
lxlO~^ Ry Bohr~^. Ultrasoft pseudopotentials together 
with a plane- wave cutoff of 35 Ry were used to describe 
the wave functions, while the charge-density cutoff was 
set to 280 Ry. The convergence threshold for the self- 
consistency of the total energy was set to 5xl0~^^ Ry, 
ensuring an accurate sampling of the complex potential 
energy surface for MOF-74-Mg. 

With the coordinates obtained from the geometry re- 
laxation, the adsorption energies and NMR shielding 
parameters were calculated using norm-conserving GI- 
PAW pseudopotentials,'^ which allow for the wavefunc- 
tion reconstruction in the atomic core region. Struc- 
tural aspects are not so sensitive to the cutoffs, en- 
suring the lower values reported above for the geom- 
etry optimizations are appropriate. But, the NMR 
shielding parameters are much more sensitive, so that 
higher values are needed. Accordingly, we used a plane- 
wave cutoff of 120 Ry and a charge-density cutoff of 
420 Ry, resulting in a convergence of the absolute shield- 
ing to within 0.05 ppm. However, in this study we 
are mostly interested in the change in shielding of the 
adsorbed molecule compared to its gas phase, which 
is converged to within less than 0.001 ppm. For the 
NMR shielding calculations,^^ we used a combination of 
the linear-response^^ and new converse^^^^ methods — 
the latter being built entirely on the theory of orbital 
magnetization.l^^H^ The adsorption energies reported in 
this study were calculated with identical parameters to 
these NMR calculations. 

In addition, we cross-checked our calculated NMR 
shieldings for selected adsorption cases of H2 with Vasf'^ 
(a plane- wave code) and Gaussian ^^ localized basis-set 
code). We find that the gas-phase shieldings of those 
codes agree to within less than 0.1 ppm with our results, 
but more importantly, the shielding difference between 
gas-phase and adsorbed molecules agrees with our calcu- 
lations to within 0.01 ppm. 



III. RESULTS 

A. Reactivity and binding energy 

Although the subject of the binding characteristics it- 
self is not the main focus of this study, we reproduce and 
extend results here that have been published elsewherd^ 
but are nonetheless important for our NMR study. In 
particular, the adsorption energy of H2, CO2, and H2O 
in MOF-74-Mg under different loading situations is of in- 
terest, as it defines the primary and secondary binding 
sites and binding geometries, for which we will report 
NMR results below. The structure of MOF-74, which 
can be seen in Figs, [l] and |2J consists of hexagonal chan- 
nels, where metal atoms at the corners are connected by 
benzenedicarboxylate linkers. The primary binding sites 



TABLE I. Adsorption energies A^ of molecules in MOF-74- 
Mg in eV for different loadings. A^^prim and AEsec are the 
average adsorption energies for the molecule at the primary 
and secondary binding sites, respectively. AE^vg is the aver- 
age adsorption energy per molecule considering all adsorbed 
molecules together. 



M 


Loading 


^-C/prim 


AEsec 


A^avg 


H2 


1 


-0.15 


n/a 


-0.15 




6 


-0.16 


n/a 


-0.16 




7 


-0.16 


-0.12 


-0.15 




12 


-0.16 


-0.12 


-0.14 


CO2 


1 


-0.50 


n/a 


-0.50 




6 


-0.50 


n/a 


-0.50 




7 


-0.50 


-0.43 


-0.49 




12 


-0.50 


-0.46 


-0.48 


H2O 


1 


-0.79 


n/a 


-0.79 




6 


-0.76 


n/a 


-0.76 




7 


-0.74 


-0.61 


-0.72 




12 


-0.74 


-0.65 


-0.70 



are located near the six metal ions in each unit cell, while 
secondary binding sites are nearer to the linkers. For fur- 
ther details, see Ref. [53l 

The adsorption energy, AE^ of a guest molecule M in 
the MOF is defined as 



AE — £^MOF+M — [^MOF + ^m] 



(1) 



where £^mof and Em are the energies of the MOF and 
the molecule in their fully relaxed form, and £^mof+m 
is the energy of the MOF with the adsorbed M. Results 
for the adsorption energies A^ are given in Table [T] for 
several different loadings: (i) low loading^ i.e. one guest 
molecule per cell occupying a primary binding site; (ii) 
high loading, six guest molecules per cell completely sat- 
urating all available primary sites; (iii) high loading, 7 
guest molecules per cell completely saturating all primary 
sites and one secondary site; and (iv) very high loading, 
12 molecules per unit cell, occupying all available pri- 
mary and secondary binding sites. For a depiction of the 
binding geometries in those cases, see Fig. [l] We find 
good agreement with the experimental adsorption ener- 
gies of -0.11 ± 0.003 eV for hJ^ and -0.49 ± 0.010 eV for 
002/^ attesting to the importance of correctly includ- 
ing van der Waals interactions in these simulations. In a 
recent study^Sl we also computed vibrational frequencies 
to obtain the thermal and zero-point energy (ZPE) cor- 
rections to these adsorption energies, allowing for a more 
accurate comparison to measured adsorption heats. But, 
for the cases considered here, we found that these correc- 
tions are on the order of 0.01 eV or less, so they are not 
again reported here. 

The influence of "crowding" on the adsorbed molecules 
in high- loading situations is present, but not dominating. 
For H2 , the contribution of the lateral interactions in the 
high-loading scenarios is negligible, being less than 7% of 




FIG. 1. The MOF-74-Mg structure is shown as sticks with metal Mg ions highUghted as green balls. The three panels show 
A) H2, B) CO2, and C) H2O as ball-and-stick representations in the case of very high loading with all six primary and all six 
secondary binding sites occupied, labelled as 1^* and 2^ , respectively. Notice the classic dimer configuration of the adsorbed 
water molecules. 



the total binding energy in the case of 12 adsorbed H2 
molecules. For the CO2 and H2O molecules, the lateral 
interactions (attractive) contribute less than 10% of the 
total binding energy as given in Table [l] when only each of 
the six primary binding sites is occupied. However, when 
all available primary and secondary sites are occupied 
(12 molecules bound), this contribution increases to 18% 
and 25% for CO2 and the hydrogen bonding H2O cases, 
respectively. 



B. NMR - Loading study 

We now move to the main topic of this study — the 
analysis of the NMR chemical shielding of H2, CO2, and 
H2O in MOF-74-Mg. Unless otherwise stated, the re- 
ported values in parts per million (ppm) are the change 
in isotropic NMR chemical shielding Aa when the gas- 
phase molecule M is adsorbed in the MOF, i. e. 



Act — (Tm in MOF — CTM in gas phase 



(2) 



After determining the primary and secondary binding 
sites for the three adsorbate molecules in the different 



loading cases (see Sec. Ill A), we calculated the NMR 
shielding of the adsorbed molecules at those positions 
within the MOF-74-Mg structure; corresponding values 
are reported in Table |TTj In primary and secondary site 
high- loading cases, the average of equivalent atoms is re- 
ported. When a single secondary site is occupied by wa- 
ter, it forms a classic water dimer by hydrogen bonding 
with the water at the nearest primary binding site. For 
this reason, only the five unpaired primary sites are av- 
eraged in Table [TI| for H2O with a loading of 7; the values 
of the "special" molecule, to which the seventh molecule 
in the secondary binding site attaches, are reported sep- 
arately. 

The chemical shielding dependence on adsorbate load- 
ing can be seen in Table [ll] for all three molecules. It 



is perhaps not surprising to see that while H2 shows a 
typically weak vdW physisorption-like interaction with 
the MOF (also indicated by the smaller binding ener- 
gies in Table |l|, water displays the more typical proton 
NMR deshielding behavior. When H2O occupies any of 
the secondary binding sites, it assumes a dimer config- 
uration with the molecule in the primary site with an 
average hydrogen-bond distance of 1.86 A which is con- 
sistent with previous results. ^^ This hydrogen bonding 
produces a large effect which can be seen by large shield- 
ing changes for H2O in high- loading cases shown in Ta- 
ble [nj These large shielding changes are in good agree- 
ment with previous NMR studies of liquid wa ter in which 
hydrogen bonding also plays a significant role .^^^^ While 
it is conceivable that as the loading increases, the water 
molecules could form a hydrogen-bond chain instead of 
occupying the remaining secondary binding sites in the 
cell, calculations and experiments^ indicate that the hy- 
drogen bond only accounts for about one third of the 
secondary site binding energy, making this configuration 
less favorable. 

From these results, it can be seen that — while direct 
usage of NMR alone to determine relative loadings for 
II2 might be difficult — the situation is more positive for 
CO2 and II2O. For CO2, relative loadings of the primary 
binding site show up clearly as a significant difference. 
While the differences are not as obvious for the secondary 
site, the presence of two peaks — one around 1 ppm and 
another greater than 2 ppm from the gas phase shield- 
ings is a good indication of CO2 occupying both primary 
and secondary sites. For II2O, the changes in shieldings 
at the primary binding site (in the absence of hydro- 
gen bonding) are much smaller as loading increases, but 
the presence of low loadings can still be detected with 
a change in shielding around 0.8 ppm less than the gas 
phase value. Furthermore, the formation of hydrogen 
bonds when greater than six molecules are adsorbed is a 
clear indication of high-loading situations for H2O. 



TABLE II. Change in NMR shielding Act in ppm (carbon 
for CO2, hydrogen for H2 and H2O) upon adsorption rela- 
tive to the gas phase for the primary and secondary binding 
sites. Shieldings for the primary and secondary binding sites 
are given in separate columns. For H2 and H2O the values 
for both hydrogens are given in separate lines. In the cases 
of H2 and H2O, the first row is the hydrogen which is more 
strongly interacting with the oxygen plane of the metal bind- 
ing site complex, while the second is further away. For higher 
loading situations, the shieldings have been averaged over all 
equivalent sites. 



^60 

CD 

E 50' 



Loading 



1 
prim. 



6 
prim. 



7 



prim. 



12 
prim. sec. 



H2 



0.48 
0.24 



0.50 
0.30 



0.54 0.09 
0.26 0.50 



0.25 0.14 
0.51 0.58 



C02 


1.01 


2.42 


2.16 


0.95 


2.69 0.81 


H20 


-0.83 


-0.85 


-i.oe^ 


' -4.41 


-0.53 -3.65 




-0.74 


-0.56 


-0.7^^ 


' -0.37 


-5.42 -0.12 



-0.57 and -4.99 ppm for the molecule to which the seventh 
H2O attaches. 




FIG. 2. Schematic drawing defining the rotation and its 
axis, depicted for the case of CO2. The molecule is placed at 
one of six identical primary binding sites, emphasized by the 
colored ball-and-stick representation with carbon, oxygen and 
magnesium represented as grey, red, and green balls, respec- 
tively. The axis of rotation is the polar axis defined as a line 
from the binding Mg through the closest atom of the small 
molecule. The molecule is then rigidly rotated (indicated by 
the curved arrow) about this axis, sampling the azimuthal an- 
gle (/? = — 360° and keeping fixed at its computed value in 
the lowest energy configuration — approximately 94° and 55° 
for H2 and CO2, respectively. 



C. NMR - Rotational study 

Strictly speaking, Table [ll] gives the NMR shielding 
change for molecules at the primary and secondary bind- 
ing sites at zero temperature. However, for finite tem- 
peratures the molecules will start to "wiggle" — governed 
by the potential energy surface around the binding site — 
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FIG. 3. The change in energy and NMR shielding calculated 
for H2 and CO2 as a function of rotation at the primary bind- 
ing site in MOF-74-Mg. For H2, the shielding for both atoms 
is shown, with the fixed hydrogen shown as solid triangles 
and the rotated hydrogen with solid circles. For a definition 
of the angle and axis of rotation, see Fig. [2] The energy scale 
is shown on the left and the NMR shielding scale is on the 
right. The dashed lines indicate ksT at room temperature, 
i. e. 25.6 meV. 



resulting in small changes in shielding. In order to inves- 
tigate the order of magnitude of such changes, we studied 
the rotational (and in the next section, translational) be- 
havior of H2 and CO2 at the primary binding site. The 
axis of rotation is defined by an imaginary line between 
the Mg atom and the closest hydrogen in H2 or the clos- 
est oxygen in CO2 at the binding site, as shown in Fig. [2] 
for CO2. The binding geometry is the lowest energy ge- 
ometry and defines the zero-degree configuration for the 
azimuthal angle (p in Figs. [2] and [3J The relative an- 
gle was then "frozen" and the molecule was rotated 
from (f = 0° to cp = 360° in intervals of 15°. Note that 
this rotation is not meant to accurately sample energies 
and shieldings of physically likely situations — rather, it 
should give an estimate for the sensitivity of these prop- 
erties in close proximity of the binding site. 

In Fig. [3] we show both the energy and NMR shielding 
as the molecules are rotated 360°. Our results show that 
for H2 at room temperature, there is a variation in the 
shielding by as much as 0.8 ppm, and at elevated temper- 
atures as much as 1 ppm. For CO2 at room temperature, 
there is a variation in the shielding up to 0.4 ppm, but at 
very high temperatures this variation can grow as large 
as 2.5 ppm. Note that there is a secondary minimum 
in the rotational energy for CO2 around 230°, so that 
low temperature measurements could detect a secondary 
peak around 2 ppm further from the gas-phase shift. This 
secondary minimum has a depth of 30 meV and thus can 
maintain trapped molecules in this rotational configura- 
tion at room temperature. 



D. NMR - Positional study 

In the previous section we investigated the sensitivity 
of energy and shielding in the proximity of the binding 
site for simple rotations. In the following, we investi- 
gate the same for large deviation from the binding site — 
in fact, we calculate the energy and shielding change 
throughout the entire cavity of the MOF. In Figs. [4J |5) 
and [6] we show a map of the energies and NMR shield- 
ings calculated for the small molecules at different posi- 
tions within the cavity of MOF-74-Mg. The values shown 
are the change in shielding when the molecules are taken 
from gas phase into the MOF in a low- loading scenario. 
The planes over which the shielding was studied are de- 
fined by three points: the center of mass of the MOF- 74- 
Mg unit cell, the binding site of the small molecule as cal- 
culated by the geometry relaxations described in Sec. [TT| 
and the Mg atom at which the molecule is adsorbed. Note 
that a plane defined in this way is not coplanar with the 
MOF structure, the normal vector to the plane being 
slightly tilted (by 10° — 15°) away from the channel di- 
rection. This implies that replicating the NMR maps 
displayed in these figures using the Dsd symmetry of the 
MOF would yield a slightly discontinuous image. The 
values were calculated at 26 equidistant points within 
this plane and linear interpolation was used to make the 
complete map. At each of the sampled points, the cen- 
ter atom of the molecule was "pinned" to the location 
of interest, and the remaining atoms allowed to relax so 
that the molecule adopted its lowest energy orientation 
and internal geometry. The row of sampled points in the 
plane closest to the MOF were found to have energies for 
all molecules of more than 2.4 eV greater than the bind- 
ing site, making it unphysical for the molecule to sample 
these regions even under high temperature and pressure 
situations. 

By definition, at the binding sites the plotted values 
coincide with the values in Table [ill Around the primary 
binding site in the energy map we have also indicated the 
25.6 meV isoline (/c^T at room temperature), which gives 
an estimate for the extent of spatial fluctuations of the 
molecules at room temperature. Note that we have also 
transferred this "fluctuation region" to the NMR maps, 
where they can now be used to estimate the fluctuations 
of the shielding change that can be expected at room 
temperature. However, it turns out that in all three cases 
these fluctuations are small. 

It is also interesting to see that even in the middle 
of the MOF all guest molecules show a shielding signif- 
icantly different from their gas phase. Although at that 
point the molecules are far away from the inside wall 
of the MOF, they are effectively surrounded by twelve 
benzene-like linkers (see the setup in Fig.fTl), the tt clouds 
of which influence the electronic structure of the adsorbed 
molecules. In the case of H2 we investigated this behavior 
further with simple model calculations using Gaussian 
with a 6-311++G(2d,2p) basis set: "Approximating" the 
MOF with 12 benzene rings and putting the H2 in the 




FIG. 4. Each triangle shows a map of the change in i) en- 
ergy and ii) the NMR shielding as the H2 molecule is moved 
throughout the cavity of MOF-74-Mg. The MOF is oriented 
so that the viewer is looking along the direction of the chan- 
nel, and hydrogen, carbon, oxygen, and magnesium atoms are 
represented as white, grey, red, and green balls, respectively. 
The two maps were calculated at identical points, but they 
are shown here in symmetry-equivalent locations within the 
MOF structure for better comparison. The shielding is aver- 
aged over both hydrogen atoms. The energy is plotted rela- 
tive to the binding energy. The white circle indicates ksT at 
room temperature, i.e. 25.6 meV, and is thus an estimate for 
the region around the primary binding site accessible through 
translational fluctuations at room temperature. 




FIG. 5. Same as Fig. |4] but here for CO2. 



center we find a difference in shielding to the gas phase 
of 0.91 ppm, whereas the full calculation with the pe- 
riodic MOF structure using QauntumEspresso gives 
0.92 ppm. 

For II2 and H2O, Figs. [4] and |6] suggest qualitatively 
similar behavior of the shielding for displacements near 
the binding site. However, as indicated in Table [ll] and 




FIG. 6. Same as Fig. ^ but here for H2O. 



Figs. [4] and [6J while H2O binds closer to the MOF (com- 
pare indicated positions of the binding sites in the maps) 
and deshields relative to the gas phase, H2 binds slightly 
further away, although when forced closer to the metal 
ion site it also deshields. 

Note that there is not a monotonic increase for either 
the energy or shielding when moving the small molecule 
from the center towards the MOF. The true function of 
shielding and energy throughout the cavity of the MOF 
is quite complicated, but we provide here a simple model 
to approximate both. To this end, we mimic the energy 
in the plane by a simple two-dimensional function (of r 
and (/) in polar coordinates) inspired by a Lennard- Jones 
potential with angular dependence, i. e. 



(ro — r)^ (ro — r)^/^ 



AE= —-- Tw^ l-bsm^f4>) +Eo. 

(3) 

Here, r is the distance from the center of the MOF and (/) 
is the polar angle. The zero angle is half-way between two 
Mg atoms. It turns out that a very similar model can also 
accommodate the spatial change in NMR shielding fairly 
well with only a minor modification, although there is no 
intuitive reason to believe the shielding should behave 
radially in a Lennard- Jones manner. For the shielding 
model, we reuse the energy model with the addition of a 
radial dependence in the angular piece of the function. 



Aa 



(ro - r)^ (ro - r) 



x/2 



l-rbsin^ {/(/)) 



do 



(4) 

These models were tested against the positional data 
used for the interpolation shown in Figs. |4]-[6J but only 
the originally sampled points were used in the parame- 
ter fitting. The parameters found for all three adsorbed 
molecules using both models are shown in Table [lllj along 
with the resulting mean absolute error (MAE) in either 
units of eV or ppm, as appropriate. 



TABLE III. Parameters for the models described in Eqs. ([3| 
and Q, along with the resulting mean absolute error (MAE), 
in either units of eV or ppm as appropriate, when compared 
to the original data. 





a 


/ 


X 


ro 


b 


Eo, ctq 


MAE 


A^ H2 


0.64 


2.79 


3.43 


6.83 


0.67 


0.15 


0.02 


CO2 


1.85 


4.84 


2.70 


6.50 


1.07 


0.48 


0.19 


H2O 


5.59 


3.12 


10.1 


6.74 


2.18 


0.86 


0.15 


Act H2 


-1.76 


-2.02 


4.24 


7.27 


0.93 


0.98 


0.19 


CO2 


10.1 


0.61 


17.8 


6.75 


-1.84 


1.01 


0.21 


H2O 


-0.19 


-1.91 


5.86 


7.03 


16.3 


1.00 


0.17 



In both models, there appear to be six free parame- 
ters, but some further simplifications can be made. First, 
there is an overall shift parameter denoted by £^0 and ao 
in Eq. ([3| and Eq. Q, respectively. These values can be 
used in the fit minimization, but essentially turn out to be 
the value of the energy change AE or change in shield- 
ing Act at the r = point in the middle of the MOF. 
Hence, these values can be fixed accordingly in order to 
reduce the number of tuneable parameters in the model. 
Second, the ro parameter in all cases is within 1 A of 
the distance from the middle of the MOF to the metal 
ion at the binding site and so could also be fixed, leaving 
only four overall adjustable parameters. We report in Ta- 
ble [Hi] these two parameters in addition to the remaining 
four so that the reader can get a sense of the magnitude 
and variance of these potentially "fixable" pieces of the 
model and more easily understand their roles in Eqs. (|3| 
and Q. 

We also find that the Aa model parameters for CO2 
reinforce what is shown in Fig. |5J i.e. that the shielding 
for CO2 adsorbed in MOF-74-Mg behaves qualitatively 
differently from H2 and H2O. There is also a difference in 
the local adsorption geometry. Whereas both the H2 and 
H2O molecules adsorb "flat-on," that is, almost parallel 
to the metal-oxygen plane which makes up the primary 
binding site, CO2 adsorbs tilted upwards as explained 
in Sec. |IIIC| pointing towards the Mg ion and interact- 
ing to a greater degree with the metal ion itself than 
with the oxygen plane. As pointed out in Ref. [79l the 
CO2 does interact with the oxygens at the binding site 
so that its rotation is affected, but it is not "pulled over" 
to the extent that H2 and H2O are. Both the quali- 
tative NMR behavior and adsorption geometry reflect a 
substantially different MOF-adsorbate interaction for the 
CO2 case from the H2 and H2O cases. 



IV. CONCLUSIONS 

In summary, we have performed ah initio DFT simu- 
lations of the energy and NMR chemical shielding of H2 , 
CO2, and H2O in the MOF-74-Mg structure. Our calcu- 
lations show that there is a relationship between loading 



in MOF-74-Mg and the change in NMR shielding. While 
the loading dependence of the shielding is small, it still 
is within the measurable experimental range. We thus 
argue that combining the NMR signal strength with the 
peak positions can yield an accurate tool for determining 
the loading of MOFs. We have further shown how the en- 
ergy and shielding behave as the molecules rotate or leave 
the binding site by providing detailed energies and NMR 
shieldings throughout the cavity of the MOF. From our 
calculated data, we were able to approximate the ener- 
gies and shieldings with two simple functions. Although 
our study only investigated one particular MOF, we be- 
lieve that the same approach is suitable for other MOFs 
and similar studies on e.g. MOF-5 and Zn(bdc)(ted)o.5 
are already in progress. 
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